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Abstract 

Background: The gene zinc finger homeobox 3 {ZFHX3) encodes a transcription factor with cardiac expression and its 
genetic variants are associated with atrial fibrillation (AF). We aimed to explore the associations between single nucleotide 
polymorphisms (SNPs) of ZFHX3 and the risk of AF in a Chinese Han population. 

Methods: 'We genotyped eight SNPs, including seven potentially functional SNPs and one previously reported SNP by using 
the middle-throughput iPLEX Sequenom MassARRAY platform. Odds ratios (ORs) and 95% confidence intervals (CIs) were 
calculated in logistic regression models. 

Results:\Ne enrolled a total of 1,593 Chinese Han origin individuals in the study, including 597 AF patients and 996 non-AF 
controls. Logistic regression analyses revealed that potentially functional SNPs rs6499600 and rsl 6971436 were associated 
with a decreased risk of AF (adjusted OR =0.73, 95% CI: 0.63-0.86, P= 1.07x1 0" 4 ; adjusted OR =0.74, 95% CI: 0.56-0.98, 
P = 0.039, respectively). In addition, rs2106261 showed a robust association with an increased risk of AF (adjusted OR = 1.71, 
95% CI: 1.46-2.00, P= 1 .85 x10 -11 ). After multiple comparisons, rs16971436 conferred a borderline significant association 
with the risk of AF. Stratification analysis indicated that the risks of AF were statistically different among subgroups of age 
for rs2106261, and the effect for rs16971436 was more evident in subgroups of patients with coronary artery disease. 

Conclusion: In summary, our study investigated the role of genetic variants of ZFHX3 in AF and two SNPs (rs21 06261, 
rs6499600) showed significant associations while rsl 6971436 conferred a borderline significant association with AF risk in 
Chinese Han populations. However, further large and functional studies are warranted to confirm our findings. 
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Introduction 

Atrial fibrillation (AF) is the most prevalent sustained arrhyth- 
mia in clinical practice characterized by uncoordinated atrial 
activation. The prevalence of AF was estimated at approximately 
2% of the general population [1,2] and this arrhythmia accounts 
for nearly one-third of hospitalizations for cardiac rhythm 
disturbances [3] . It is closely related to a 5-fold risk of stroke, a 
3-fold incidence of congestive heart failure, and higher all-cause 
mortality, which may contribute to high healthcare cost (average 
2000 EUR per AF patient per year) [4-6]. 

AF frequently occurred coupled with multiple cardiovascular 
risk factors such as advancing age, male sex, hypertension, 
diabetes, ischemic, valvular heart disease, heart failure, and 
hyperthyroidism [7-9]. However, some patients maintain sinus 
rhythm regardless of various risk factors for AF while some other 
patients develop AF without any evidence of AF-related diseases 
[10]. Epidemiologic studies have shown that first-degree relatives 



of AF patients were 1.77-fold to 4.67-fold more likely to have AF 
than the general population, and offspring with one parent 
suffered from AF had approximately a 2-fold increase in the risk of 
developing AF [11,12]. These facts indicate that predisposing 
genetic factors may affect AF susceptibility. Over the past years, 
lots of genetic studies of AF were performed by applying several 
techniques such as linkage analysis, candidate gene resequencing 
and association studies. Using these methods, especially genome 
wide association studies (GWASs), investigators have identified 
multiple genes and genetic loci which are associated with AF 
[13,14]. In 2009, Benjamin et al. [15] conducted a meta-analysis 
of GWASs for AF and identified a locus for AF (ZFHX3, 
rs2106261, RR = 1.19; P= 2.7 xlO" 7 ). Candidate gene studies are 
approaches employed by researchers to screen variants in any 
gene that is hypothesized to cause a disease. We sought to identify 
potentially functional AF susceptibility SNPs of £FHX3 using the 
candidate gene tool. Furthermore, we aim to confirm the 
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Table 1. Clinical characteristics in AF cases and AF-free controls. 



Variants 


Cases (N = 597) 


Controls (N = 996) 


P value 


Male gender (%) 


397(66.5%) 


674(67.7%) 


0.630 


Age, years 


58.4±11.5 


59.0±10.2 


0.278 


Paroxysmal AF (%) 


383 (64.2%) 


NA 




Persistent AF (%) 


196 (32.8%) 


NA 




Permanent AF (%) 


18 (3.0%) 


NA 




Lone AF (%) 


71 (11.9%) 


NA 




Hypertension (%) 


260 (43.6%) 


267 (26.8%) 


<0.001 


Diabetes (%) 


53 (8.9%) 


28 (2.8%) 


<0.001 


CAD (%) 


48 (8.0%) 


51 (5.1%) 


0.019 



Data are presented as mean ± standard deviation or number (percentage). 
AF, atrial fibrillation; CAD, coronary artery disease; NA, not available. 
doi:1 0.1 371 /journal.pone.01 01 31 8.t001 



association of rs2 106261 in ^FHX3 with AF risk in Chinese Han 
populations. 

Materials and Methods 

Study population 

This case-control study was approved by the Ethical Committee 
Review Board of Nanjing Medical University, China. All the 
participants enrolled in the study gave their written informed 
consent. 

In brief, the cases were incident AF patients recruited from 
Department of Cardiology, the First Affiliated Hospital of Nanjing 
Medical University from June 2010 to August 2013. All AF cases 
in the study were confirmed according to the diagnostic criteria, 
2011 ACCF/AHA/HRS focused updates incorporated into the 
ACC/AHA/ESC 2006 guidelines for the management of patients 
with atrial fibrillation [3] . Physical examinations were carried out 
by expert cardiologists using routine 12-lead electrocardiography 
(ECG), and/or ambulatory ECG recordings. According to clinical 
characteristics, AF can be classified into paroxysmal AF (episodes 
that generally last 7 days or less), persistent AF (episodes that 
sustain beyond 7 days) and permanent AF (ongoing long-term 
episodes, in which cardioversion has failed or has not been 
attempted). Young AF individuals (under 60 years) without clinical 
or echocardiographic evidence of cardiopulmonary disease, 
including hypertension were considered as lone AF [16]. The 
control subjects were genetically unrelated inpatient individuals 
from multiple departments of the First Affiliated Hospital of 
Nanjing Medical University. All the controls were confirmed to be 
free of AF based on ECG or medical files at the time of 
enrollment. 

We also collected general and clinical information from medical 
recording files in the hospital system such as age, gender, and 
history of hypertension, diabetes, coronary artery disease (CAD), 
and hyperthyroidism. Patients with hyperthyroidism, severe 
cardiac dysfunction (NYHA Class IV) and advanced age (beyond 
90 years) were excluded in both AF and control groups. All the 
enrolled individuals are of the ethnic Chinese Han origin by self- 
report. 

SNP selection 

We first used public HapMap SNP database (phase II + III Feb 
09, on NCBI B36 assembly, dbSNP bl26) to search SNPs that 
localized within the gene region of %FHX3 (including 10 kb up- 



stream region of the gene), with MAF>0.05 in Chinese Han 
population (CHB). Then, a web-based analysis tool was used to 
predict the function of these SNPs (http:/ /snpinfo.niehs. nih.gov/ 
snpinfo/snpfunc.htm). After function prediction analysis, a total of 
9 potentially functional SNPs were selected. Using the Haploview 
4.2 software, linkage disequilibrium (LD) analysis with an r z £0.8 
was further applied to filter these functional SNPs and 7 
(rsl2596992, rsl3336412, rsl6971312, rsl6971436, rsl7680796, 
rs6499600, rs8049936) of which were remained. Additionally, 
previously reported single nucleotide polymorphism (SNP) 
rs2 106261 in £FHX3 was selected as well for this study. As a 
result, 8 loci were finally determined to perform genotyping. 

SNP genotyping 

Blood samples were drawn from study participants and genomic 
DNA was isolated from EDTA-preserved whole blood, using the 
standard phenol-chloroform method [17]. The genotyping was 
performed by using the middle-throughput iPLEX Sequenom 
MassARRAY platform (Sequenom, Inc, San Diego, CA, USA). 
All SNPs were successfully genotyped with call rates >95%. 

Statistical analysis 

The standard independent samples i-test was used to compare 
continuous variables and the % 2 test was used to compare 
categorical variables. The Hardy-Weinberg equilibrium compar- 
ing the observed genotype frequencies with the expected ones was 
determined by the % 2 test in the control group. Odds ratios (ORs) 
and 95% confidence intervals (CIs) were calculated using logistic 
regression analyses in the additive model to assess the strength of 
the associations between the variants genotypes and AF risk. The 
heterogeneity of associations between subgroups was assessed 
using the x 2 -based Qjtest. All statistical analyses were performed 
using STATA 12.0 software (Stata Corporation, College Station, 
TX, USA) and P<0.05 was considered significant. 

Results 

Characteristics of the Study Population 

A total of 1,593 Chinese Han origin individuals were recruited 
in this study, including 597 AF patients and 996 non-AF controls. 
The demographic characteristics of cases and controls are 
summarized in Table 1. As shown, there were no significant 
differences for the distributions of age and gender between cases 
and controls (P values were 0.278 and 0.630, respectively). The 
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mean age of AF and control groups was 58.4±11.5 and 
59. 0± 10.2 years, respectively. Male proportion was 66.5% and 
67.7% in AF group and control group, separately. Among cases, 
64.2% patients were paroxysmal AF, 32.8% were persistent AF, 
while others were permanent AF. And patients with lone AF 
accounted for a small part (11.9%) of all AF cases. Considered as 
common risk factors for AF, hypertension, diabetes, and CAD 
were more prevalent in cases than controls (P<0.05). 

Overall Associations between ZFHX3 Variants and AF risk 

The genotype distributions of the 8 polymorphisms and their 
associations with AF were shown in Table 2. The observed 
genotype frequencies for these SNPs were all in agreement with 
Hardy-Weinberg in the controls. Logistic regression analyses 
revealed that functional SNPs rs6499600 and rs 1697 1436 were 
associated with decreased risks of AF in the additive model 
(adjusted OR = 0.73, 95% CI: 0.63-0.86, P= 1.07 x 10" 4 ; adjusted 
OR =0.74, 95% CI: 0.56-0.98, P = 0.039, respectively). In 
addition, rs2 106261 showed a robust association with an increased 
risk of AF (adjusted OR =1.71, 95% CI: 1.46-2.00, 
P = 1.85x10 ). We further calculated Pvalues for false discovery 
rate (P-FDR) to perform multiple comparisons. After comparisons, 
we found that rs6499600 and rs2 106261 remained their associ- 
ations with AF risk while rsl6971436 lost its association with the 
risk of AF. There were no obvious evidences of significant 
association between other 5 SNPs and AF risk. 

Stratification Analysis 

In the stratification analysis, we further evaluated the associa- 
tions of rs6499600, rs 1697 1436 and rs2 106261 on AF risk in 
subgroups based on age, gender, hypertension, diabetes and CAD. 
As shown in Table 3, a significant difference among subgroups of 
age was observed for the association of rs2 106261 with AF risk (P 
for heterogeneity =0.001). Patients with CAD exhibited more 
significant protective effect than patients who were free of CAD for 
rsl6971436 (adjusted OR =0.13, 95% CI: 0.02-0.71). 

Discussion 

In this study, we investigated the associations between 
polymorphisms of ZTHX3 and the risk of AF in a case-control 
study including 597 cases and 996 controls. Seven potentially 
functional SNPs and one previously reported polymorphism in 
ZFHX3 were selected for genotyping in our population. Two 
SNPs (rs6499600 and rs2106261) were identified to be significandy 
associated with AF risk and rsl6971436 conferred a borderline 
significant association with the risk of AF in Chinese Han 
populations. Forrs2 106261, the effects were statistically different 
among subgroups of age in stratified analysis, and the risk 
associated with rsl6971436 was more evident in patients with 
CAD. 

^FHX3 (also known as ATBF1), encodes a transcription factor 
with multiple homeodomains and zinc finger motifs, and regulates 
myogenic and neuronal differentiation. Little is known about the 
function of %FHX3 for cardiac expression, but it has been shown 
to interact with Protein Inhibitor of Activated Stat3 (PIAS3) [18]. 
STAT3 was found to be a regulator of paracrine circuits in the 
heart, and it was essential for interstitial matrix deposition balance 
and capillary vasculature maintenance [19]. Increased expression 
of STAT3 has been observed in animal models of AF and 
proposed to contribute to atrial matrix deposition. Tsai et al. 
suggested that activated angiotensin II/Racl/STAT may be 
associated with or perhaps contribute to the structural and 
inflammatory changes in AF [20] . 
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A total of 3 SNPs were significantly associated with the risk of 
AF (rs2 106261, rs6499600 and rsl6971436) in our population. 
The SNP rs2 106261 locates in the first intron of ^FHX3 and has 
been reported to confer increased risk of AF in European and 
Chinese populations [15,21]. Even though the mechanism of 
rs2 106261 conferred increased susceptibility of AF is little known, 
one plausible hypothesis is that it may alter the expression of 
ZTHX3 by modulating regulation of ^FHX3 transcription or pre- 
mRNA splicing because intronic sequences have been found to 
play regulatory roles in gene expression recently [22]. The SNP 
rs6499600 is located in the seventh intron of £FHX3. And it may 
affect the binding of transcription factors such as ATF6 (activating 
transcription factor 6), which was predicted by the online tool 
SNPinfo (http:/ /snpinfo. niehs.nih.gov/index.html). ATF6 is a key 
transcriptional activator to maintain cellular homeostasis, and 
inhibition of which induced dilatation of left ventricle and 
depression of cardiac function in murine heart [23]. ZTHX3 
rsl6971436 is a T to G variant at the codon of 428 of exon 1, 
resulting in an amino acid alteration from threonine (Thy) to 
proline (Pro). Therefore, we proposed that the missense variant 
might lead to the change of 2J<HX3 expression, consequendy 
cause a dysregulation of STAT3 and finally result in altered 
susceptibility of AF. 

Several limitations of the present study need to be addressed. 
Firstly, we did not replicate the results in additional individuals, 
this may contribute to potential false positive errors. The present 
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analysis was restricted to individuals of Chinese Han descent, and 
therefore, the findings may not be generalizable to individuals of 
other races and ethnicities. Secondly, as echocardiography results 
were unable to be obtained from every participant, left ventricular 
ejection fraction (LVEF) values were not measured to evaluate 
cardiac function in both groups. Also, we did not analyze the 
impact of SNPs on stroke because some patients enrolled in the 
study had not performed head CT or MRI scans. Furthermore, 
the way to select potentially functional SNPs as target ones by 
using web-based tool may bring some positive or negative errors. 
Therefore, the results are required to be further replicated by well- 
designed studies in additional large-scale Chinese Han popula- 
tions. 

In conclusion, our study investigated the role of genetic variants 
of ZFHX3 in AF and two SNPs (rs2 106261, rs6499600) showed 
significant associations while rs 1697 1436 conferred a borderline 
significant association with AF risk in Chinese Han populations. 
However, further large and functional studies are warranted to 
confirm our findings. 
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